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On the Presence of Thermal SZ Induced Signal in the First 
Year WMAP Temperature Maps 
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ABSTRACT 

Using available optical and X-ray catalogues of clusters and superclusters of galax- 
ies, we build templates of tSZ emission as they should be detected by the WMAP 
experiment. We compute the cross-correlation of our templates with WMAP temper- 
ature maps, and interpret our results separately for clusters and for superclusters of 
galaxies. For cluster s of g alaxies, we claim 2-5 a detections in our tem plates built from 
BCS jlEbeling et alJT99al200Ch . NORAS llBohringer et alJ200nh and lde Grandi et alJ 
( 1999j) catalogues. In these templates, the typical cluster temperature decrements in 
WMAP maps are around 15-35 fiK in the RJ range (no beam deconvolution applied). 
Several tests probing the possible influence of foregrounds in our analyses demonstrate 
that our results are robust against galactic contamination. On supercluster scales, we 
detect a diffuse component in t he V fc W WMAP bands w hich cannot be generated by 
superclusters in our catalogues ijEinasto et alJll99"3.l!997j) . and which is not present in 
the clean map of Tegmark, de Oliveira-Costa & Hamilton (2003) . Using this clean map, 
our analyses yield, for Einasto's supercluster catalogues, the following upper limit for 
the comptonization parameter associated to supercluster scales: ysc < 2.18 x 10~ 8 at 
the 95% confidence limit. 

Key words: cosmic microwave background - galaxies: clusters: general - diffuse 
radiation - intergalactic medium 



1 INTRODUCTION 

The WMAP experiment is currently observing the mi- 
crowave sky at five diff erent frequency chann els with differ- 
ent angular resolution feennett et al]l2003al) . The low fre- 
quency bands, K, Ka & Q, measuring at 23, 33 & 41GHz 
respectively, are expected to be particularly sensitive to the 
free-free and synchrotron emission in the Milky Way. Al- 
though their angular resolution is not as good as in the 
high frequency channels, their measurements of the fore- 
ground contamination are critical in order to achieve an 
optimal subtraction of non-cosmological signal in the over- 
all analysis. This analysis has yielded the identification of 
two main components in the sky signal: a cosmological 
component, which constrains crucial cosmological param- 
eters such as Ha, S7 m , w, T re j orl with unprecedented accu- 
racy iBennett et aljl2003blh and a foreground-induced com- 
ponent, whose impact in the high frequency channels (V 
(61GHz) & W(94GHz)) is modelled by means of the low fre- 
quency measurements. These channels have been located in 
a frequency range where the contribution from foregrounds 
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is expected to be minimal jBennett et alJl2003dh and their 
high angular resolution (up to 0.21°) enables the study the 
sub-horizon structure of the Last Scattering Surface. 

However, it is also expected that the presence of Large 
Scale Structure intersecting the geodesies of the CMB pho- 
tons leave a signature in the form of secondary temperature 
anisotropics. Among these temperature fluctuations, the 
most important (in terms of amplitude) are those caused 
by the Integrated Sachs- Wolfe effect, (ISW), due to the 
variation of the gravitational potent ial, and the thermal 
Sunya ev-Zel'dovich (hereafter tSZ, ISunvaev fc Zeldovicbl 
( 1980)), due to the distortion of the black body spectrum of 
CMB photons after interacting with a hot electron plasma. 
The search of tSZ-induced signal in CMB data sets by com- 
paring them with surveys of Large Scale Structure has been 



perfo rmed for C OBE data (Gan ga. Cheng. Mever. fc Page 
19931: iBouehn fc JahodaT 11991 Bandav et al 



199S iKneissl et alJ Il997lh and Tenerife data 
^^bm&^arthi^^djio^Tjarandela. fc Hernandez-Monteagudol 
2000). They all failed to measure any significant correlation, 
and hence could only set upper limits to the comptonization 
parameter y. 

In this context, three recent analyses have been car- 
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ried out on WMAP data. This first one has been carried 
out by the WMAP team: a cross-correlation of the W- 
band with the XBA Cs catalog of 242 Abell-type clusters 
febeling et al Il996l) has prompted a detection at 2.5a level. 
However, a direct cross-correlation of WMAP data with 
ROSAT ijDiego. Silk fc Sliwal Etioll l has given negative re- 
sults. In this analysis, the difference map between V and 
W bands was compared to the ROSAT map by the study of 
the power spectrum of the product map and the cross-power 
spectrum of both maps. Apart from an apparently spurious 
coincidence in the product map, th ere is no trace for any 
signifi cant correlation. More recently. iBoughn fc Crittenden! 
(2003) have reported a positive cross-correlation at the 2-3 
sigma level between the WMAP data and the X-ray HEAO- 
1 and radio NRAO VLA sky surveys. Their study focuses 
on extra-galactic objects as tracers of the Large Scale Struc- 
ture, and the positive sign of their correlation would point 
to the ISW effect as the r esponsible for this excess. A s im- 
ilar result was reported by Fosalba & Gaztanagal ll2003l) af- 
ter cross-correlating the V band of WMAP with a template 
built from the APM galaxy survey. They find a positive cross 
correlation with a significance of ~ 2a. 

In this work, we extend the study of possible correlation 
of WMAP data with Large Scale Structure (LSS) to cluster 
and supercluster scales. By using cluster surveys present in 
the literature, we construct template maps of LSS, and cross 
correlate them with the temperature maps of WMAP. Al- 
though the angular resolution and the sensitivity of WMAP 
are not ideal for the typical tSZ amplitudes and angular pro- 
files of clusters of galaxies, we expect that if a large enough 
number of them are indeed contributing to the map their fi- 
nal imprint should be statistically measurable. On the other 
hand, angular sizes of superclusters of galaxies are resolved 
by the V & W bands of WMAP, so they should allow to set 
constraints in the density and/or temperature of the diffuse 
gas thought to be present in those structures. 

In Section 2 we outline the details of the statisti- 
cal methods used in our analyses. Section 3 describes the 
WMAP data, whereas in Section 4 we explain how our tSZ 
templates of clusters and supercluster of galaxies were syn- 
thesised. Finally, in Section 5 we show our results, which we 
discuss in Section 6. 



2 STATISTICAL METHOD 

In this section we describe the statistical tools used to 
search for signature of Large Scale Structure in WMAP 
CMB maps in the form of tSZ effect. All our statistics 
will be defined in the real space, and will be restricted to 
those regions where we expect a high contribution from SZ 
sources. By doing this, we are trying to minimise the effects 
related to the limited angular resolution and sensitivity of 
WMAP in terms of typical SZ amplitudes. A Fourier analy- 
sis in a selected and non-continuous sample of pixels cannot 
always be defined in a simple way and will be discarded 
here. Core radii of clusters of galaxies usually subtend an 
angle of ~ l' — 3' in the sky, which is considerably smaller 
than the FWHM of the highest resolution channel (94GHz) 
of WMAP, (~ 12.6 arcmins). For this reason, from the 
five bands of WMAP, we shall only use this channel when 
searching for cluster induced signal in the WMAP data, as 



we expect that other channels will dilute the SZ signal to a 
much greater extent. The high resolution map of Tegmark, 
de Oliveira-Costa & Hamilton (2003), which combines all 
bands but preserves the angular resolution of the 94GHz 
channel, will also be considered. However, for superclusters 
of galaxies, the situation differs considerably, provided the 
fact that these sources can have a size of a few degrees on the 
sky, enabling the use of lower frequency channels of WMAP. 

The two statistical methods used in this paper will 
search for similarities present in both our templates and the 
WMAP CMB maps. The first method (method i) consists 
of a pixel-to-pixel comparison of CMB maps to the tem- 
plates. Given the CMB and noise amplitudes in those pixels, 
we estimate what is the fraction of signal present in CMB 
maps which is correlated to the structure of our catalogues. 
In the second method, (method it) we computed the cross- 
correlation function between our templates and the WMAP 
map(s), and compared it to the auto-correlation function 
computed in our synthesised maps. 

We shall assume that the WMAP data can be modelled 
as the following sum of components: 



T = T cmb + aM + ^ ^F 1 + N, 



(1) 



where T C mb denotes the CMB, N the noise, and M is the 
SZ induced signal modelled by our templates, which enters 
into the map with an average amplitude given by a. CiF* 
denotes the ith foreground component. This last term will 
introduce a bias in our statistics, and its impact in our 
results must be considered. There might be as well some 
offset due to the inaccurate sett ing of the zero level in the 
WMAP maps. From figure (7) in lBennett et alJ (I2003J) one 
would expect a typical error of a few (~ 3 — 5) microkelvin. 
The impact of such offsets will be addressed when quoting 
our results. 



Method i): 

We shall minimise the statistic x 2 given by 



X 



(T-aM)C _1 (T-aM) 5 



(2) 



in terms of a. By doing this, our estimate of a will be equal 
to: 



E[a] 



TC M 



(3) 



MC-!M T 

In both equations|5|and|21 C denotes the correlation function 
of the SZ-free temperature map, i.e.: 



C —< (T — aM) 1 (T - aM) > 



(4) 



In our analyses, we considered only CMB and noise when 
computing C. The cosmological signal was modelled by using 
the CMB power spec trum measured by the WMAP team 
(Ben nett et alJ l2003aT) . whereas the noise component was 
introduced using the amplitudes at each pixel according to 
the scanning strategy of the WMAP mission 1 . The formal 
error for the estimated a is given by: 



1 All data related specifically to the WMAP mission has been 
obtained from: |http: / /lambda.gsfc.nasa.gov| 
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However, in these equations we have neglected the bias 
term introduced by the presence of any possible foreground 
component, (i.e. synchrotron, free-free, dust or point 
sources). Our basic assumption will be that foregrounds 
and tSZ sources will introduce temperature fluctuations 
of opposite sign, positive and negative respectively, in the 
WMAP CMB scans. We remark at this point that all 
WMAP channels are in the Rayleigh- Jeans (RJ) frequency 
range, for which the tSZ distortions of the CMB planckian 
spectrum introduce negative temperature fluctuations. 
Therefore, because of the presence of foreground residuals 
in WMAP maps, the measured a will be shifted to positive 
values, and we will only be able to establish lower limits to 
the presence of SZ-induced signal in WMAP scans, so our 
approach can be regarded as conservative. 

Method ii): 

We shall define the \ 2 statistic in terms of the cross- 
correlation and auto-correlation functions computed from 
the WMAP CMB map(s) and our templates. In this case, 
we write \ 2 m the form: 



\- =^[Ctm(^) 



■ aCMM(0i)] 



^[CTM^-aCMMie,)] (6) 

In this equation, Ctm and Cmm denote the cross and 
auto correlation function of our template maps with the 
WMAP scans, respectively. 7r is the covariance matrix, de- 
fined as: 



< (Ctm(0 4 ) - <C TM (0i) ) ) (Ctm(^) - {c TM (e 3 ) ) ) ). (7) 

The angle brackets denote ensemble averages after com- 
puting the cross-correlation between our templates and 100 
Monte-Carlo realizations of the microwave sky accounting 
for the CMB signal and instrumental noise. The value of a 
which minimises \ 2 an d its formal error are given by the 
following equations: 



E[a] 



J2ij C ™(8i)' n ' i j Cmm{0j) 
C M M{Qi)^C M M{Qj , 



y ' CMM{6i)l^ i j'CMM{0j 



(8) 



(9) 



This method assumes that all components which are 
not correlated to the template M average to zero when 
computing the cross-correlation. Foregrounds do not aver- 
age to zero, and so they are supposed to bias the estima- 
tion of a towards positive values. Therefore, the impact of 
foregrounds in this analysis is quite similar to that in the 
previous method. We must also note that any mismatch in 
the zero level of the map would yield spurious correlation, 
so it is convenient to make alternative tests, such as repeat- 
ing analyses after rotating the templates, to check whether 
any measured cross correlation is actually associated to our 
templates or not. 

We expect method ii to be particularly sensitive to 



the spatial structure of superclusters of galaxies. Indeed, 
under WMAP's 94 GHz beam, almost all clusters can be 
regarded as point-like, so we should expect Cmm{0) = 
for 9 > Ofwhm- However, superclusters of galaxies can 
subtend up to 40 degrees in the sky, and for them it might 
be possible to trace their typical angular extent by the study 
of Ctm (0). For these reasons, we shall use method i for 
the search of tSZ signal in cluster scales, whereas method 
ii will be applied on supercluster templates. Due to the 
typical size of superclusters (few degrees), one can include 
in the analysis different bands of WMAP, as we show in the 
following section. Nevertheless, it is worth to remark that, in 
terms of the frequency dependence of the tSZ effect, the sig- 
nal should drop around a 14% from channel V to channel W. 



3 THE WMAP DATA 

We have mentioned that, due to angular resolution require- 
ments, initially we considered only the W band for our anal- 
yses with templates of clusters of galaxies. Nevertheless, for 
the sake of comparison, we included the map provided by 
Tegmark, de Oliveira-Costa & Hamilton (2003), here de- 
noted as T03. This map weights the five bands of WMAP 
according to their sensitivity at every angular scale. The re- 
sult is a map with as much angular resolution as the W band 
and smaller presence of both noise and foreground contam- 
inants. However, the properties of its noise component are 
complex. In our analyses, we characterised it in a optimistic 
and/or in a conservative way: either we assumed that the 
excess power of the map (er excess — 66/iK) was due to gaus- 
sian white noise 2 , or we simply assigned WMAP's 94GHz 
band noise pattern to T03. 

In order to assess the impact of foregrounds on our re- 
sults, we performed our analyses using two different fore- 
ground masks: the KpO mask and the Kp2 mask, both pro- 
vided by the WMAP team. These masks remove from the 
analysis those pixels showing temperature fluctuations that 
trespass a given threshold computed from the histogram 
of the K-band data. The digit in the name of the mask 
is proportional to the threshold, so KpO is twice as con- 
servative/strict as Kp2. In these analyses, the HEALPix 3 
resolution parameter was N s ide =512. 

The relatively large typical size of superclusters allowed 
us to test the consistency of our results by including the V 
band in our analyses. Furthermore, clean maps (in the sense 
of foreground-free) provided by the WMAP team (here 
denoted as ilc for internal linear combination map) and by 
Tegmark, de Oliveira-Costa & Hamilton (2003) were also 
considered. The ilc map has a typical angular resolution of 
one degree, whereas T03 has as much angular resolution 
as the W map. However, at low angular resolution, both 
maps are very similar, so we only retained T03 in our 
analyses. For our purposes of studying superclusters, the 
resolution of one degree suffices, so those temperature maps 
and all supercluster templates where pixelized under a 
coarser resolution, of a pixel size of ~ 56'. This angular 

2 This gaussian noise was assumed to follow the spatial pattern 
of the noise templates provided by the WMAP team 

3 HEALPix URL site: http://www.eso.org/science/healpix/ 
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degradation allowed us to speed the computation of the 
cross-correlation between the WMAP temperature maps 
and our supercluster templates, and also made feasible 
to estimate the covariance matrices of equation (J7J under 
a reasonable amount of time. We expect that the linear 
combination of different bands used to build the T03 map 
does not remove the tSZ component of our maps, but leaves 
it in roughly the same amplitude level. Our results will 
justify this hypothesis. 



4 TEMPLATES 

One crucial step in our procedure will be constructing the 
templates from existing surveys of clusters and superclus- 
ters of galaxies, both in the optical and in the X-range 
band. We shall build only tSZ templates, that is, in units 
of decrements of thermodynamic temperature. When cross- 
correlating with WMAP temperature maps, a positive cor- 
relation would point to a detection of tSZ effect. In any case, 
what matters is the spatial structure of our templates and 
not their overall normalisation. The conversion into thermo- 
dynamic temperature decrements from X-ray based cata- 
logues will be performed by using existing relations between 
flux and/or luminosity in the X-ray band and antenna tem- 
perature decrements. For the optical catalogues, further as- 
sumptions will be necessary. 

In the case of galaxy clusters, there is a fairly large 
sample of catalogues available, both in the X-range and in 
the optical. In the northern sky, our analyses will first fo- 
cus on two X-ray based catalogues of clusters of galaxies, 
i.e., the No rthern ROSAT All Sky Galaxy (NORAS) Clus- 
ter Survey feohringer et al.ll200ri) and the ROSAT Bright- 
est Cluster Sample (BCS. lEbeling et alJ lll998t |2000|U . 
The source of these two catalogues is t he publish e d dat a 
of the ROSAT All Sky Survey (RASS, iTriimperl il993ft . 
IVogues et all jTjjj) ). The NORAS cluster survey contains 
495 sources with extended X-ray emission, of which 378 are 
unambiguously identified as clusters of galaxies. Due to the 
modest completeness of this catalogue, (count-rate limit of 
0.06 counts/s), we extend our analysis to the BCS catalogue, 
which is flux limited and shows a completeness of around 
90%. This catalo gue is built on two di fferent cluster samples: 
the first sample (lEbeling et aljH998Tl l isted 201 clusters o f 
galaxies, whereas the extended sample (lEbeling et al]l2000T) 
provided 99 new objects. The criteria used in their construc- 
tion are not purely X-ray based. Both catalogues (NORAS 
and BCS) give luminosities in the 0.1-2.4 keV energy band 
for every cluster, which will be converted into brightness 
temperature fluctuations, as it will be shown. 

In the southern hemisphere, Ide Grandi et alJ <ll999l) 
extracted from the RASS data a flux limited sample of 
130 clusters of galaxies, with an estimated completeness of 
around 90 %. Redshifts and fluxes in the 0.5-2 keV energy 
band are given for every object. 

Analysing R OSAT Position Sensit ive Proportional 
Counter pointings, IVikhlinin et alJ <ll998T) compiled a cat- 
alogue of 203 serendipitously detected galaxy clusters. This 
catalogue, hereafter denoted by V98, covers both hemi- 
spheres, and therefore complements all catalogues listed 
above. V98 provides redshifts and fluxes of members in the 



Table 1. SZ cluster sample used to calibrate the template maps 



Cluster Name 


z 


AT™ M L {0-l-2AkeV) (b) 
(mK) (10 44 h~Q erg s" 1 ) 


A1656(COMA) 


0.0232 


-0.55 ± 0.10 


7.26 


A2256 


0.0601 


-0.44 ±0.09 


7.11 


A2142 


0.0899 


-0.90 ±0.14 


20.74 


A478 


0.09 


-0.92 ± 0.15 


13.19 


A1413 


0.143 


-0.96 ± 0.11 


13.28 


A2204 


0.152 


-0.96 ± 0.28 


21.25 


A2218 


0.171 


-0.75 ± 0.20 


9.30 


A665 


0.182 


-0.91 ± 0.09 


16.33 


A773 


0.197 


-0.89 ± 0.10 


13.08 


A1835 


0.252 


-1.34 ±0.15 


38.53 


Z3146 


0.291 


-0.86 ± 0.14 


26.47 



(a) Central decrement, from the compilation of lCoora^ ^ 

(b) X ray luminosities from the BCS catalogue. 



energy range of 0.5 - 2 keV. A conversion into luminosities 
will be necessary in order to estimate the antenna temper- 
ature fluctuations induced by these clusters. All catalogues 
listed so far rely on ROSAT data, and hence are all condi- 
tioned by the observation strategy of this experiment. 

For the sake of comparison, we shall include in the anal- 
ysis two more catalogues of galaxy clusters. These remain- 
ing two catalogues were built on the ba sis of existing opti- 
cal galaxy cluster surveys, namely APM (jPalton et alll997f) 
and ACO (lAbell. Corwin fc Qlowinll989l l. These catalogues 
provide an homogeneous sample of galaxy clusters in both 
hemispheres, and have been processed after applying selec- 
tion criteria well different from those used in X-ray data 
analysis. However, unlike in the previous cases, these cata- 
logues are not directly sensitive to the hot gas causing SZ 
distortions, and assumptions will have to be made when re- 
lating optical properties of the sources (i.e. richness) with 
antenna temperature decrements. Given the angular resolu- 
tion of the W band, we note that we can safely ignore issues 
related to the gas distribution within cluster scales. 

With respect to superclusters, we shall recu r to the cat- 
alogues provided bv lEinasto et alJ (Il994l Il997ll . The use of 
these templates will be essential when testing the hypothesis 
of the presence of diffuse hot gas distributed in megaparsec 
scales. 

The last step in our template construction procedure 
was assigning temperature decrements to the sources in our 
catalogues. For the case of galaxy clusters in X-ray base d 
catalogues, it was done in the same way than lCooravNl999lh 
by relating AT t sz with Lx from a sample of clusters with 
measured RJ antenna temperature decrements 4 Following 
his discussion, we would expect this relation to have the 
form ATtsz = aL x . A \ 2 minimisation 5 was used to derive 



4 It should be noted that quoted values correspond to the in- 
ferred central temperature. These numbers are generally derived 
by fitting a /3-profile to t he observational data (see individual ref- 
erences at lCoora^<1999D l. so the telescope dependence is reduced, 
and a model dependence is introduced. 

5 We must remark th at we are using here a different statistic 
than in lCooravl ll999l) . so our error bars are different. 
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the best-fit pair (a,b), using luminosities in the 0.1-2.4 keV 
ROSAT band from the BCS catalogue 

/ A0.l-2.4ksV) \0^±0.09 

AT t % = -(0.24±0.06) —f— 2 r mK{10) 

\ 10 ho er 9 s 1 J 

The values used for this fit are summarised in table Jy. 
A similar fit can be obtained from the luminosities quoted 
in the NORAS catalogue (a = 0.21 ± 0.05, b = 0.53 ± 
0.09). Clusters extracted from the RASS-PSC catalogue had 
fluxes in the same energy range than NORAS and BCS 
(0.1 - 2.4 keV), so a conversion into flux and ATtsz was 
straightforward. For V 98 and the southern catalogue of 
Ide Grandi et alJ (I1999T) . a previous scaling of fluxes from 
the 0.5 - 2 keV energy band into the 0.1 - 2.4 keV band 
was first required (we obtained an average conversion fac- 
tor of Fa (0.5 - 2keV)/F x (0.1 - 2AkeV) » 0.65). For op- 
tically selected clusters (ACO and APM catalogues), we 
have assumed that the observed SZ decrement is propor- 
tional to some power n of the cluster richness R. Follow- 
ing feennett et al]|l9 93h. we have probed here their values 
n — 1,2,3. In addition, we have used the cluster sample 
described in table Q to calibrate a relation with the form 
L x = A(R + 1)™' , obtaining 



- (0.1-2.4feeV) 



(1.6±0.3)(R+1) 1A±0 - 2 10 44 hso erg s _1 (ll) 



We finally related Lx with ATtsz by using the scalings given 
above. This template will be labelled as ACO calibrated. 

We conducted two different approaches when modelling 
the templates of superclusters of galaxies. We first intro- 
duced an uniform distribution of hot gas centred on the 
optical center of the supercluster and with radius the an- 
gular size of the supercluster, as it is provided by the cata- 
logues. This template will be referred to as Sph-SC for Spher- 
ical Superclusters, and will probe the presence of hot gas in 
SC scales. However, some correlation present between these 
structures and the WMAP CMB maps might be diluted due 
to the huge angular size (9sc ~ 40°) of some members. For 
this reason, our second approach consisted in filling with hot 
gas a circle of radius one degree around every cluster mem- 
ber in each supercluster. In some occasions, these circles 
overlap along a given direction in those superclusters where 
member clusters align forming a filamentary structure. We 
shall label this template as rldg-SC. We assumed that the 
spheres of hot gas in both Sph-SC and rldg-SC templates 
had the same density and temperature in every case, and so 
we assigned identical tSZ decrements to all of them. 

There were two types of sources in our templates, 
according to their angular size. Point-like sources refer to 
all those clusters which are not resolved by the W channel 
of WMAP 6 . Extended sources refer to COMA, VIRGO 
and superclusters. In all cases, the templates have been 
convolved with a realistic approximation of the average 
WMAP 94GHz beam. In this initial stage, the pixelization 
of all templates was such that the typical pixel size was 
around 7', (N si d e = 512 under HEALPix pixelization). 
In the case of superclusters, the templates were degraded 
down to N a ide = 64, which corresponds to a pixel size 
of ~ 54'. Therefore, for superclusters, templates do not 



Some non-resolved clusters are seen in the W map, as A2142. 



contain 3, 145, 728 pixels, but just 49, 152 pixels. 



5 RESULTS 

5.1 Cluster templates 

In the present section we search for correlations between 
our cluster templates and WMAP CMB scans. After build- 
ing the cluster templates, they were convolved with the sym- 
metrised window functions corresponding to each of the four 
differencing assemblies (hereafter DA's) present in the W- 
band. These convolved templates were weighted according 
to the noise pattern of each DA, yielding a final template 
which was produced in exactly the same way as the averaged 
CMB map corresponding to the W-band of WMAP. In order 
to investigate which clusters contributed to the CMB maps, 
we applied three different amplitude masks, so only pixels 
brighter than some threshold would be used when searching 
for correlation. In practice, it consisted in computing the 
histograms of pixel amplitudes for each template. These his- 
tograms typically showed a symmetric distribution around 
zero amplitude (due to the numerical noise introduced by 
the convolution) plus a bright tail, formed by cluster pixels. 
By increasing/lowering the threshold amplitude we merely 
change the number of pixels in the tail considered in the 
analyses. For the three masks of increasing amplitude, here- 
after labelled as z/i, v% & 1/3, we selected a sample of 2000, 
500 and 50 pixels, respectively. 

We used four different CMB maps in this section. One 
of them was the W band of WMAP. This channel con- 
tains a white noise of an average amplitude of 175/iK, which 
makes it significantly more noisy than the map provided by 
Tegmark, de Oliveira-Costa & Hamilton (2003). For com- 
parison, we included this map (T03) in our analyses as well. 
Furthermore, in order to check for systematic effects in our 
approach, we repeated the analyses after rotating both maps 
5 degrees around the axis orthogonal to the Galactic plane, 
i.e., 51 — 5°. The rotated maps will be labelled as R-W 
and R-T03 for W and T03, respectively. By comparing re- 
sults from W and T03 to R-W and R-T03, we can i) check 
how our error bars compare with the case of lack of correla- 
tion, and ii), check that any detected correlation vanishes in 
the rotated map, so that we can assign such correlation to 
sources of angular size < 5° . As mentioned earlier, we try to 
track the impact of foregrounds by performing our analyses 
in both KpO and Kp2 masks. 

All results are shown in table for the KpO mask. In 
this table we define (5 as /3 = a(M), where M denotes the 
template. That is, j3 is given in temperature units (fiK), and 
contains the sign of the cross-correlation, i.e., negative for 
tSZ. From this table, one can deduce that three catalogues 
are giving some significant tSZ induced signal, namely BCS, 
NORAS and de Grandi. For all of them, the statistical sig- 
nificance of the detection is higher under the column T03 7 . 



7 Actually, no detection is reported for deGrandi in the W band. 
Nevertheless, if one takes into account the error associated to 
noise in this band, it turns out that the values of (3 under W are 
compatible (always below the 2a level) to those quoted in the 
T03 column. 
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This is partially associated to the fact that in this case, we 
are taking the optimistic values of the noise level for the T03 
map. The pessimistic approach would consist in quoting the 
values of j3 under the T03 column but with the error bars 
given in the W column. We must remark, however, that the 
error bars quoted in the optimistic case are not far from (but 
rather well within) the null (i.e., free of tSZ emission) values 
of /3 quoted under the R-T03 column. With respect to op- 
tical catalogues, we can appreciate that the APM catalogue 
at n = 2, 3 in T03 gives values for j3 which are consistently 
increasing with increasing v, surpassing the threshold of 2a 
in n = 3, v\. This trend might be clarified with the second 
year WMAP data release, for which noise amplitude is ex- 
pected to be a factor of v2 lower. In all cases where we claim 
detection, ATtsz increases with the threshold, although its 
statistical significance diminishes as a consequence of the 
smaller amount of pixels. The only slight increase of tSZ 
amplitude with the threshold assures that the correlation is 
raised not exclusively by the very brightest clusters, but by 
most of the clusters considered in each catalogue. Otherwise, 
the signal would be much more diluted with decreasing v s . 
Note that these amplitudes give the average contribution of 
the tSZ effect to WMAP maps, and have to be corrected 
for the beam and pixel smearing in order to give the actual 
cluster decrements. These analyses were repeated with zero- 
mean CMB maps, yielding similar results. The presence of 
possible offsets in the CMB maps therefore has negligible 
impact on our conclusions. 

5.2 Supercluster templates 

We next describe the results obtained after applying our 
method ii on the V, W & T03 maps. We restricted our 
analyses of correlation functions to a maximum angle of 
dmax = 80°, as this is roughly twice the largest size of the 
sources in our catalogues. In figure Q we show the cross 
correlation functions for KpO and rldg-SC. Shaded regions 
are given by the Monte Carlo simulations and show the un- 
certainty region around (TM) associated to instrumental 
noise and cosmic variance. For the T03 case, we used the 
optimistic approach to characterise the noise field. From the 
similarity of all shaded areas in our three bands, we conclude 
that instrumental noise cannot be of importance at the an- 
gular scales probed by our supercluster template. We can 
see that for the V & W bands, the rldg-SC template gives 
a positive correlation with the CMB scans. Both bands give 
roughly the same amplitude for the cross-correlation, and 
are also very close to those obtained under the Kp2 mask. 
However, if the template is rotated 45 degrees in the di- 
rection of galactic longitude, the cross-correlation does not 
drop, but remains in roughly similar levels. Testing the pos- 
sibility of a simple offset in the zero-level between the V & 
W bands and T03, this analysis was repeated after subtract- 
ing the mean of the Supercluster templates. In this case, the 
cross correlation function was compatible with zero for all 
angles. Moreover, this result is also obtained after applying 
our method ii on the T03 map: the cross-correlation func- 
tion is not significantly different from zero in the 6 range 

8 If only pixels at V3 are responsible for the correlation, then one 
would expect /3 to drop by a factor of 10, 40 for U2, ux, respectively. 



of [0°,80°] 9 We know that the T03 map is supposed to be 
cleaned from foreground signal to a much greater extent 
than the V & W bands. And we have also seen in our anal- 
ysis of clusters that the tSZ signal was not removed when 
combining the different frequencies to build the T03 map, 
at least in the cluster scales. For the large scales (low multi- 
poles) we have also checked from figure (6) of Tegmark, de 
Oliveira-Costa & Hamilton (2003) that the tSZ signal is pre- 
served. Therefore, we associate the correlation present in the 
V & W bands to either an offset or a diffuse foreground com- 
ponent present in these bands but mostly absent in the T03 
map. Accordingly, we shall try to give physical interpreta- 
tion to those results obtained from the T03 map exclusively. 

In the T03 map, the Sph-SC catalogue gave very sim- 
ilar results to those of rldg-SC, (see figure ©)• Our sta- 
tistical analysis (method ii) failed to report any significant 
detection in any case, but allowed us to limit the comp- 
tonization parameter associated to these angular scales to 
bey < 1.86 x 10" 8 ,2.18 x 10" 8 for Sph-SC under the KpO, 
Kp2 masks, respectively, at 95% c.l. The rldg-SC template 
yielded y < 8.59 x 10 -9 for both masks, and again at 95% c.l. 
However, one must keep in mind that our modelling of the 
Supercluster structure has been rather simplistic, and that 
these limits do not observe the inaccuracies of our model 
when describing the density and temperature distribution 
of the diffuse gas. 

5.3 Ring analysis 

In previous subsections, we have demonstrated that there is 
evidence of tSZ- induced signal detection using X-ray cluster 
templates, while we do not have such detection either for 
optical cluster templates or for supercluster catalogues. In 
order to check the robustness of these numbers, we have also 
perform ed a "ring" analysis, simi lar to that uti l ized p revi- 
ously in lBennett et al.l i 199311 and lBandav et all Jl99fj) . We 
then compute the difference between the temperature in the 
pixel where we expect to have a cluster, and the mean tem- 
perature in a ring of radius between 9i and 62- We used 
here 6»i =1.5FWHM, and 6 2 = v^i- With this analysis, we 
should reproduce similar amplitudes for the above quoted 
detections, with the advantage that now we are particularly 
unsensitive to a possible zero offset of data. The obtained 
values using the KpO mask and the W & T03 maps are shown 
in table Error bars take into account the temperature 
correlation of the CMB component among different pixels. 
From this table, we confirm the detection of BCS, NORAS 
and de Grandi (this one only in the T03 map) catalogues. 
The obtained numbers are comparable to those values of 13 
using V2 (which roughly corresponds to use the central pixel 
of all clusters in each catalogue), and the signal dissappears 
when analysing the rotated maps (R-W and R-T03). The 
other catalogues (V98, Vogues, ACO and APM) provide no 
detection, although we marginally see the aforementioned 
trend in APM using T03. 

In order to perform a similar test to our study on su- 
percluster scales, we have applied a ring analysis in which 

9 At 9 ~ 60° we found an excess of anticorrelation which dissa- 
peared after removing the quadrupole in the CMB temperature 
map. 
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(3±a p , QuK), [KpO] 



Template Catalogue 


v\ 


w 

20.34 ± 6.83 


R-W 
18.57 ± 6.84 


T03 
-1.32 ± 4.75 


R-T03 
5.21 ± 4.82 


ACO (L x =L X {R+1)) 


V2 


17.21 ± 11.26 


13.23 ± 11.32 


-2.40 ± 7.60 


2.72 ± 7.64 




v-i 


32.35 ± 28.69 


12.46 ± 28.54 


24.36 ± 18.10 


-1.49 ± 18.27 




VI 


13.44 ± 6.65 


13.89 ± 6.64 


0.54 ± 3.78 


-0.35 ± 3.84 


ACO (n=l) 


VI 


16.22 ± 10.80 


11.31 ± 10.85 


-3.41 ± 7.03 


-0.08 ± 7.14 




V3 


20.69 ± 28.13 


-23.49 ± 27.96 


5.36 ± 17.29 


-35.93 ± 17.56 




V\ 


6.93 ± 6.02 


9.83 ± 6.05 


-0.54 ± 3.13 


-2.72 ± 3.22 


ACO (n=2) 


VI 


2.76 ± 10.91 


4.72 ± 11.07 


-6.81 ± 6.27 


-3.99 ± 6.54 




V3 


16.47 ± 28.41 


-29.00 ± 28.50 


5.19 ± 17.36 


-30.98 ± 18.00 




V\ 


-0.07 ± 4.79 


4.36 ± 4.91 


-2.40 ± 2.26 


-2.26 ± 2.41 


ACO (n=3) 


VI 


0.66 ± 10.23 


2.34 ± 10.44 


-3.81 ± 5.53 


-3.33 ± 5.91 




V3 


-24.07 ± 30.24 


-34.18 ± 30.92 


-22.77 ± 17.74 


-10.30 ± 19.30 




V\ 


-3.37 ± 8.94 


-2.22 ± 8.93 


-5.67 ± 4.16 


-3.89 ± 4.15 


APM (n=l) 


V2 


-2.19 ± 15.96 


6.15 ± 15.95 


-7.57 ± 10.40 


-11.41 ± 10.42 




V3 


3.77 ± 31.99 


29.28 ± 32.05 


-3.05 ± 21.92 


2.82 ± 22.09 




V\ 


-5.58 ± 5.74 


1.57 ± 5.76 


-3.24 ± 2.55 


-3.69 ± 2.58 


APM (n=2) 


VI 


-10.63 ± 12.45 


9.96 ± 12.45 


-9.62 ± 6.27 


-9.10 ± 6.30 




V3 


-0.54 ± 31.73 


36.01 ± 31.78 


-13.10 ± 19.49 


-3.33 ± 19.95 




V\ 


-4.16 ± 3.39 


-0.01 ± 3.43 


-3.52 ± 1.48 


-1.50 ± 1.52 


APM (n=3) 


VI 


-6.62 ± 9.07 


5.25 ± 9.17 


-5.93 ± 4.10 


-7.19 ± 4.22 




V3 


-20.71 ± 30.95 


30.18 ± 31.09 


-17.91 ± 18.20 


-8.31 ± 18.84 




V\ 


-3.75 ± 3.05 


-0.32 ± 2.98 


-7.24 ± 1.32 


-0.14 ± 1.29 


de Grandi 


VI 


2.10 ± 10.45 


2.75 ± 10.27 


-21.49 ± 5.19 


5.31 ± 5.01 




V3 


6.37 ± 28.68 


-27.48 ± 28.47 


-48.58 ± 17.72 


25.55 ± 17.41 




v\ 


-18.10 ± 5.30 


6.61 ± 5.37 


-8.44 ± 2.50 


-2.76 ± 2.57 


BCS 


V2 


-29.87 ± 11.40 


16.66 ± 11.62 


-19.00 ± 6.63 


-8.25 ± 6.81 




V3 


-64.06 ± 31.78 


19.59 ± 31.96 


-39.30 ± 18.43 


-12.30 ± 18.62 




Vl 


-10.95 ± 5.47 


8.12 ± 5.44 


-10.47 ± 2.64 


4.12 ± 2.69 


NORAS 


V2 


-15.04 ± 11.06 


18.17 ± 10.97 


-17.89 ± 6.50 


3.79 ± 6.60 




v-i 


-29.65 ± 29.16 


42.09 ± 28.51 


-49.54 ± 17.30 


1.50 ± 17.37 




Vl 


5.52 ± 4.31 


2.98 ± 3.92 


0.24 ± 1.89 


-0.06 ± 1.73 


V98 


V2 


12.94 ± 12.20 


11.11 ± 11.22 


-0.56 ± 6.68 


1.18 ± 6.11 




V3 


17.87 ± 33.01 


26.42 ± 29.55 


-1.32 ± 19.89 


-5.13 ± 17.75 




Vl 


0.50 ± 2.44 


1.12 ± 2.37 


0.34 ± 1.07 


0.19 ± 1.05 


Vogues 


V2 


14.91 ± 8.78 


11.10 ± 8.61 


3.47 ± 4.07 


0.86 ± 4.05 




V3 


45.18 ± 31.03 


37.81 ± 30.22 


10.09 ± 18.51 


1.57 ± 18.23 



Table 2. This table shows the results of method i on cluster templates for the mask KpO for the real beam. /3 is defined from a 
as /3 = c?(M), it is given in fiK and contains the sign of the cross-correlation, (i.e., /3 < in the presence of tSZ induced signal). 
Similar results are found for mask Kp2. Each template has been cross-correlated to the W band of WMAP, to the combined map of 
WMAP produced by Tegmark, de Oliveira-Costa & Hamilton (2003), T03, and to both maps rotated by 5° in galactic longitude, R-W 
& R-T03, (see columns). For each template, we show in different rows results for different thresholds applied on the templates. More 
severe threshold (more luminous clusters) give higher amplitudes for the mean AT t gz m those catalogues where a significant detection 
is reported, namely, BCS, NORAS and de Grandi, although the latter shows significant detection only in the T03 map. Note that the 
optimistic estimates error bars for T03 are not too far from the null detections listed in the R-T03 column. 



we substract the mean temperature in annuli of inner/outer 
radii [8i =1.5FWHM, 62 = 1°] from the mean temperature 
in annuli limited by the inner/outer radii [62, O3 = \/202], in 
all pixels where we have a cluster member of the superclus- 
ter. As for our meth od ii, we find no trace of t SZ emission 
due to diffuse gas in lEinasto et al.l il994L IT997I) catalogues, 
nor in cluster or superclusters scales. 

Before ending this subsection, we would like to remark 



the agreement between methods i and ii and the ring anal- 
ysis. 



6 DISCUSSION AND CONCLUSIONS 

We have used optical (Einasto, APM, ACO) and X-ray 
based (BCS, NORAS, de Grandi, V98, Vogues) catalogues 
to locate in the sky tSZ emitting sources as clusters and su- 
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Table 3. Results of the ring analysis on cluster templates for 
the W and T03 maps, using the KpO mask. We also include the 
results derived rotating the maps 5° (R-W and R-T03). All values 
are in fiK. 



KpO 



Catalog 


w 


R-W 


T03 


R-T03 


BCS 


-28.1 ± 11.0 


-5.6 ±10.9 


-15.3 ± 5.0 


±0.3 ±4.9 


NORAS 


-16.6 ±8.6 


±4.5 ± 8.6 


-9.9 ±3.9 


-0.7± 3.9 


Grandi 


-2.4 ± 16.2 


-1.0 ± 16.0 


-27.3 ± 7.4 


±4.5 ± 7.3 


V98 


+4.6 ± 12.6 


+2.2 ±12.5 


-0.9 ±5.7 


±4.6 ± 5.7 


Vogues 


-6.3 ±13.3 


-7.4 ±13.3 


-11.2 ± 6.1 


-7.7 ±6.1 


ACO 


±4.3 ±2.7 


±5.2 ± 2.7 


-0.3 ±1.2 


+0.3 + 1.2 


APM 


-5.1 ± 5.6 


±3.0 ±5.6 


-3.5 ±2.6 


-0.2 + 2.6 


KpO 



t 4 




Figure 1. This figure shows the cross correlation function (TM) 
computed from our rldg-SC template and the V (thick solid line), 
W (dotted line) and T03 (thick dashed line) CMB maps. Each 
line is centred in a shaded region limiting the l-cr confidence re- 
gion computed from Monte Carlo realizations. We associate the 
high amplitude of the cross correlation for the V & W bands 
to foregrounds. T03 gives results compatible with zero. The auto 
correlation function {MM) amplified by a factor of 10 is displayed 
by the thin solid line. Analyses under the Kp2 mask yielded very 
similar results. 



perclusters of galaxies. Using scaling relations present in the 
literature, we have assigned microwave fluxes to our sources, 
and built tSZ catalogues as they would be seen by the W 
band of WMAP. We have cross correlated these catalogues 
with WMAP using two different statistical approaches for 
clusters and superclusters. For the BCS, NORAS and de 
Grandi cluster catalogues, we claim detections at the 2-5 a 
level, (except for deGrandi, which gives no detection in the 
W band). In these cases, we have checked that most of the 
clusters (and not only the brightest) are contributing to the 
correlation, with amplitudes of typically 20-30 fjK in the 
WMAP scans. We fail to detect any statistically significant 
signal from the (optical) ACO and APM catalogues, (al- 
though the latter has chances to give a detection if the noise 





9 (degrees) 

Figure 2. Same as in Fig.lTt. but for the Sph-SC template. 



level decreases after the second year of observations). This 
may reflect the fact that optical catalogues are not sensitive 
to the temperature of the IGM in clusters, and might include 
many clusters that, being present in the optical, contribute 
with negligible distortions to the tSZ signal, hence diluting 
the overall cross-correlation. We believe the impact of fore- 
grounds to be of negligible relevance in our results, provided 
the fact that our method yields very similar results for dif- 
ferent masks (KpO and Kp2) and different CMB scans, (W 
band of WMAP and the clean map of Tegmark, de Oliveira- 
Costa & Hamilton (2003, T03)). Moreover, all these results 
are recovered when a ring analysis (consisting in studying 
the temperature fluctuations in annuli centered on cluster 
pixels) is applied on our catalogues. 

With respect to our supercluster templates, we have 
tested the hypothesis of hot diffuse gas comptonizing the 
CMB photons in megaparsec scales. We have built two differ- 
ent supercluster templates, depending upon we place the hot 
gas uniformly in the superclusters or concentrated around 
the cluster members. By computing the cross correlation 
function between these templates and the CMB scans, and 
comparing it with the template auto-correlation function, 
we have traced the level at which supercluster induced sig- 
nal is present in the CMB maps. We find that results de- 
pend largely on the CMB map we cross-correlate with: the 

V & W WMAP bands give positive cross-correlation, several 
cr-levels above zero. This correlated signal does not vanish 
when the template is rotated 45 deg in galactic longitude, 
but drops to zero in the T03 map. The T03 map keeps the 
same zero temperature level of all bands of WMAP, but 
combines them in a multipole dependent manner to min- 
imize the presence of foregrounds. Hence, we point to an 
offset or foregrounds as the most likely responsibles for the 
different cross-correlation present between our supercluster 
templates and the W & T03 CMB maps, and we discard the 

V & W bands from our cross-correlation function analyses. 
When applying our method ii on T03, we find no significant 
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cross-correlation , and this allows us to place the follow- 
ing constraint on the supercluster induced comptonization 
parameter ysc associated to the Einasto supercluster cata- 
logues: ysc < 2.18 x 10 -8 at the 2a or 95% confidence level. 
These figures, however, do not account for the uncertainty 
related to the modelling of the gas distribution. Finally, we 
conclude noting that any use of the V & W bands in cross- 
correlation analyses is at best compromised by the presence 
of a diffuse foreground component, at least at the level of a 
few microkelvin. 
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" The ring analysis ap plied on supercluster catalogues of 
lEinasto et alj Jl994, 1997) also fails to detect any tSZ induced 
signal. 



